Efficient fold-change detection based on protein-protein interactions 
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Various biological sensory systems exhibit a response to the relative change of the stimulus, often 
reffered to as fold-change detection. Here, we present a mechanism consisting of two interacting 
proteins, able to detect a fold-change effectively. This mechanism, in contrast to other proposed 
mechanisms, does not consume chemical energy and is not subject to transcriptional and transla- 
tional noise. We show by analytical and numerical calculations that the mechanism can have a fast, 
precise and efficient response for parameters that are relevant to eukaryotic cells. 

PACS numbers: 



I. INTRODUCTION 

A fold-change detecting (FCD) mechanism with input 
u and output y is characterized by the property 



At 



(1) 



indicating that the output signal depends only on a rel- 
ative change in the input signal rather than the absolute 
value FCD fits various biologial sensory systems, such 
as vision and hearing. Recent studies show that FCD 
mechanisms are present in E. coli and Xenopus laevis 
embryos ||. 

In the last few years, it was demonstrated that the in- 
coherent feedforward loop [1,(3 and the two-state protein 
Q can provide FCD. As a consequence of the involvement 
of continues protein production and degradation, they 
consume chemical energy, and are subject to transcrip- 
tional and translational noise 0- Besides, transcription 
and translation set a lower limit for the response time, 
which can be quite high, especially in mammalian cells 
(~ 20min) @-[l0j]. 

Here, we describe a FCD mechanism, based purely 
on protein-protein and protein-ligand interactions. The 
mechanism it is not subject to transcriptional and trans- 
lational noise and does not consume chemical energy 
since it does not involve protein production and degra- 
dation. The characteristic timescales are set by the rates 
of protein interactions, which have a much broader range 
than the transcription and translation rates fllT - fli} . 

This paper contains three sections. Section [TT1 describes 
the mechanism and shows formally that the mechanism 
can provide FCD. In section Hill the action of the mech- 
anism is demonstrated by numerical simulation. In this 
section, it is shown that the mechanism can be preciese, 
robust and efficient for biologically relevant parameters. 
Section IIVI provides a discussion and outlook of the re- 
sults. 



II. THE MECHANISM 

The mechanism we suggest uses two proteins, named 
X and Y. The system responds with the output Ys to 
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FIG. 1: A schematic representation of the FCD mechanism 
as proposed in this paper. The mechanism has input S and 
output Ys- For clarity, Y is drawn two times. The copies are 
grouped by a box. 



the concentration of the input ligand, named S 1 , where 
Ys is the complex of Y and S. The mechanism involves 
three reactions, as presented schematically in Fig [TJ In 
the first reaction, 5* and X bind, forming the complex 
Xg. The binding and unbinding rates are given by fc™ s 



and , respectively: 
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In the second reaction, S and Y bind, forming the output 
complex Ys with binding rate ky^ and unbinding rate 

toff. 
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In the third reaction, X$ and Y bind, forming the com- 
plex XgY with binding and unbinding rates given by 
fc™ s y and k^ sY , respectively: 
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X s + Y <=* X S Y. 
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(4) 



Fig. [T] shows a schematic representation of the mecha- 
nism. 



A. The mechanism satisfies FCD conditions 

In general, a mechanism based on input S, internal 
variable X, and output Ys, described by functions 



X s = f(X s ,Y s ,S), 
Y s =g(Xs,Y s ,S). 



provides FCD if [|: 



f(pX s ,Ys, P S)=pf(Xs,Y s ,S), 
g(pX S7 Y,pS)=g(X s ,Y s ,S), 



(5) 
(6) 



(7) 
(8) 



for any p > 0. Below, we show that the suggested mecha- 
nism can provide FCD. This is accomplished by showing 
that conditions (0 and (j8]) can be satisfied. 
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The first condition, 0, requires the ratio S/Xg for an 
input $, which increases by a factor 0, 



5(i) = 



5o if t < to 
4> So if t > to 



(9) 



to be independent of the initial input concentration Sq. 
When concentration Y is small compared to concentra- 
tions S and X , the major reaction that involves X$ is the 
binding of X and S 1 given by reaction © . The mean- field 
concentration X$ as a function of time is described by the 
Michaelis-Menten equation fT" 
can be approximated by 



which for k°^ X < kf 



X, 



k%XS 



kf s X s . 
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By comparison with condition ([7]), it follows that, in this 
case, the first condition is satisfied. 



FIG. 2: Mean-field concentrations for S = lOnM, X = fiM, 



Y = lOraM, = 10'M _ V 



fcx" = 10" 



b on — 



lOPM -1 *) -1 , kf s = 10s" 1 , fc<f g y = lO^M -1 *- 1 , and kf sY - 
lOs" 1 as a function of time. Here, InM corresponds to 1 
protein for prokaryotes and 10 3 proteins for eukaryotes. The 
input S fold-changes at t = 1.5s and t = 4s by a factor 4> — 3. 



the Michaelis-Menten equation that the equilibrium value 
of Ys is approximated by 



Ya = 



S/k Ys 

Xs/kxsY 



Y. 



(11) 



where k XsY = kf sY /k™ Y and k Ys = kfjk™. By 
comparison with condition ©, it follows that, in this 
case, the second condition is also satisfied. 



2. Second condition, ((8} 

The second condition, ©, requires the output Ys to 
be proportional to the ratio S/Xs- When reactions ([3]) 
and (01 equilibrate fast compared to reaction @, a con- 
tinues equilibrium for reactions involving Y can be as- 
sumed. If Y has a strong affinity to Xs, it follows from 



B. Mean-field FCD 

In sum, it was shown that the mechanism provides 
FCD if 

a. concentrations S and X are large compared to Y; 

b. reaction rate fc™ X is small compared to k 1 ^ ; 
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FIG. 3: Mean-field solution and random walk simulation of protein concentrations for S — IQnM, X — fiM, Y — WriM, 



kT„ = 10 7 A/- 
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= 0.1s" 



fc$? = 10 9 M" 



I.off 



= 10s"\ k°£ sY = 10 8 M _1 s _1 , and = 10s" 1 . The input S is 

presented on the upper plots. It fold-changes at t = 2, 4.5, 7s and t = 9.5s by a factor = 3. On the bottom plots mean-field 
solution is represented by the thick gray lines while thin black lines represent random walk simulations. For eukaryotes, the 
two lines mostly overlap. Here, InM corresponds to 1 protein for prokaryotes and 10 3 proteins for eukaryotes. 



c. protein Y has a strong affinity to X,s compared to S; 

d. reactions and (U]) equilibrate fast compared to re- 
action ([2]). 

When the conditions are satisfied, the FCD transient re- 
sponse is characterized by two timescales, 



n 



k^ n s YS 



and T2 



X s 



(12) 



corresponding to the raising and the decay of the re- 
sponse, respectively. 

Fig. [5] shows an example of mean field values X, Xs, 
Y, Yg, and XgY as a function of time for an input given 
by J9]) with = 3 and to = 1.5s and t = 4s. One can see 
that the response, Yg, does not depend on the absolute 
value of the input, S. This shows that themcchanism 
provides FCD for this set of parameters. The parameter 
values and the robustness to noise are discussed in Sec. 

nil 



III. EXAMPLES AND BIOLOGICAL 
RELEVANCE 

Affinity constants of protein-protein and ligand- 
receptor interactions span, over six orders of magnitude, 
1CP 12 — 10 _6 M, a result of the large range of the as- 
sociation and dissociation rates which are typically in 
the 10 3 - 10 10 M -1 s -1 and 1CT 4 - 10 4 s _1 range, respec- 
tively [IlT - [l4| ]. Protein concentrations typically vary in 



the range 1 — 10 3 nM, corresponding to 10 3 — 10 6 pro- 
teins for eukaryotes and 1 — 10 3 proteins for prokatyotes 
[l6| . Below, we show that such a large range of biochem- 
ical parameters enables fulfilment of the FCD conditions, 
still leaving a lot of freedom. This freedom, in principle, 
can be used, for instance, to tune the timecale of the 
detection response. 

In order to demonstrate the relevancy of the mecha- 
nism to biological systems, we analyzed equations ([2] [3l 
2]) using mean-field and numerical random walk calcula- 
tions [l8| . Fig. [3] shows the result for a certain set of 
parameters for prokaryotes and eukaryotes, which differ 
by the number of proteins present. This parameter set 
is used also in Fig. [5J The mean-field response peaks 
of the system are approximately transient and identical, 
which shows that the mechanism can provide FCD. Note 
that the last peak is, in this case, only approximately 
transient since condition (|b|) is no longer satisfied. For 
prokaryotic cells, the protein number is too low for an 
effective FCD on the level of a single cell. For eukary- 
otes, the output signal is nearly unaffected by the finite 
number noise. In this example, the timescales of the ini- 
tial response is less than a tenth of a second; the decay 
timescale is approximately 0.25s. 



IV. DISCUSSION AND OUTLOOK 

In this paper we propose a new FCD mechanism, 
based purely on protein-protein and ligand-protein inter- 
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actions. The mechanism, in contrast to previously pro- 
posed mechanisms [1, 0] , does not require continues tran- 
cription, translation and degradation of proteins. The 
mechanism avoids transcriptional bursts and other noise 
sources due to small number bottleneck processes Q , and 
benefits from the typical large numbers of proteins in eu- 
kariotic cells. The initial response time, in principle, is 
limited only by diffusion of proteins [llj], and, therefore, 
can be smaller than a second for a biologically relevant 
set of parameters, as shown in Fig. [3] 

For eukaryotes, the noise due to finite number statistics 
does not play a significant role. However, for prokary- 
otes, the typical protein number is not sufficient for a 
precise FCD on a single cell level. Previously proposed 
mechanisms are relatively consuming, since the number 
of proteins has to be large and their lifetime has to be 



short in order to reduce a detection noise and make the 
responce quick. From this perspective, the mechanism 
as proposed here is advantageous, since it is acting in 
detailed balance [l9| . 

We showed that the mechanism is characterized by two 
timescales, t\ and T2 given in Eq. (fT2|) . t\ characterizes 
the raise of the response while T2 characterizes its decay. 
The value of t±, in principle, is limited only by diffusion of 
proteins [12| . Thus, as shown in the example, presented 
in Fig. it can be smaller than a second for a biologically 
relevant set of parameters. 

We expect that the proposed mechanism is only a sin- 
gle example of a new, large class of biochemical fold- 
change detectors, based purely on ligand-protein and 
protein-protein interactions. 
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